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We examine experimentally the deformation of flexible, microscale helical ribbons with nanoscale 
thickness subject to viscous flow in a microfiuidic channel. Two aspects of flexible microhelices 
are quantified: the overall shape of the helix and the viscous frictional properties. The frictional 
coefficients determined by onr experiments are consistent with calculated values in the context of 
resistive force theory. Deformation of helices by viscous flow is well-described by non-linear hnite 
extensibility. Under distributed loading, the pitch distribution is non-uniform and from this, we 
identify both linear and non-linear behavior along the contour length of a single helix. Moreover, 
flexible helices are found to display reversible global to local helical transitions at high flow rate. 


Helices have captured the fascination of many for cen¬ 
turies, from Darwin’s observation of plant tendrils ^ to 
a child’s play with a toy Slinky. Beyond curiosity, the 
interaction of small helices with fluids is particularly im¬ 
portant because of its relevance to both fundamental sci¬ 
ence mg and technological applications, such as swim¬ 
ming microrobots or microflow sensors [TOHig . Nature 
has perhaps best demonstrated the importance of small 
scale helix-fluid interactions through the evolution of he¬ 
lically shaped flagella, which are exploited by swimming 
microorganisms to move through their surrounding flu¬ 
ids [MIH]. At these length scales, structures function at 
low Reynolds number (i.e. inertia is negligible and vis¬ 
cous forces play a dominant role), defined as Re = pvl/ri, 
where p and ry are the fluid density and viscosity, re¬ 
spectively, V is the flow velocity, and I is a characteristic 
length scale. In these instances, the helical structure is 
key to locomotive functionality; however many questions 
remain with regard to the fluid-helix interactions at these 
small length scales. 

While helices in low Reynolds number flows have been 
considered in several studies over the past couple decades, 
experimental work has focused mainly on macroscopic, 
non-deformable helical models in high viscosity fluids 
gnang, likely due to the difficulties in fabricating and 
analyzing microscopic systems in a controlled manner. 
A natural bacterial flagellar filament is on the order of 
tens of nanometers in diameter and several microns long 
with bending stiffness in the range of B ~ 10“^'^ to 10“^^ 
N m^ [2OH22] . values that have been measured through 
optical tweezer or crude flow experiments. This low flex¬ 
ural stiffness results in drastic changes in shape of bac¬ 
terial flagella observed experimentally under the motion 
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of fluids [ini [T^ Uni 1131 m] ■ Moreover, the frictional co¬ 
efficient that defines the relative resistance of motion be¬ 
tween the solid and fluid is an important physical param¬ 
eter for small helices in flow H [MHl Hi ng. Therefore, 
a microscopic experimental model that examines flexible 
helices in low Reynolds number flow, with the ability to 
predict and extract helical shape changes and frictional 
properties, would be exceedingly beneficial. 

In this Communication, we examine the deformation of 
synthetically fabricated helical ribbons in controlled vis¬ 
cous flow with length scales and mechanical properties 
that approach those found in bacterial flagella and mi¬ 
croscale robots [Mini ( i.e. microscale radius, nanoscale 
thickness). Here, taking the helical ribbon thickness 
~50 nm as the characteristic length scale for axial flow 
experiments, our approach allows for experiments in low 
Reynolds number even with strong flows (^10 mm/s). 
Re ^ 10“^. We discuss our findings in the framework of 
resistive force theory [MI17] and demonstrate that, as 
expected, the size, shape and bending stiffness of a heli¬ 
cal ribbon defines the axial deformation of microhelices 
in flow [ZlEej. We quantify the non-uniform shape of a 
flexible helix deformed by viscous drag, showing that the 
pitch distribution transitions from linear to non-linear 
behavior within the same helix as a function of flow ve¬ 
locity. Our measurements allow us to assess validity of 
the resistive force theory and extract the effective fric¬ 
tional coefficient for microscale, flexible helices. 

Consider a helical ribbon defined by its axial length 
{H), contour length (L), pitch (p), and radius {R), as well 
as its cross-section, which is defined by the ribbon width 
(w) and thickness (t), as illustrated in Fig. To cre¬ 
ate such structures, we recently reported a method that 
relies on spontaneous formation of helices from initially 
flat ribbons, driven by 2-phase elastocapillary deforma¬ 
tion HH]. The ribbons are taken to be inextensible {i.e. 
a fixed contour length), and under the condition that 
t/w 1 and w/L <C 1, helices form by bending in the 
direction of the nanoscale thickness (Figs. and[^). 
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FIG. 1. (a) Geometry of a helical ribbon, (b) The relationship between the thickness and width of PMMA ribbons is best fit 
to a quadratic: t = aw^ where a = 3530 m“^. Inset: 3D profile of ribbon cross-section measured by optical profilometry. y-axis 
is in nm and x-axis is in Note the nanoscale thickness and microscale width, (c) Experimental setup that allows helical 
ribbons to form in a large pool and be placed into a connected microfluidic channel and (d) the placement of the helix in the 
vertical center of the channel (at Vmax)- (e) Measured stretch ratio and flow rate in a 3 cycle experiment with corresponding 
fluorescent images in (f). 


A key point to emphasize is that the preferred helical 
radius has a strong dependence on the ribbon thickness 
[5S]; hence the bending stiffness, B — El and the helix 
radius, R are not independently controlled {E being the 
Young’s modulus and I ^ wt^ being the second moment 
of area). This approach is advantageous since it provides 
versatility in controlling the helix geometry through con¬ 
trol of fabrication parameters. 

In our experiments, a flow rate (Q) is applied to a he¬ 
lix that is held in a microfluidic channel. To fabricate 
the helices, ribbons are first prepared on a flat substrate 
by an evaporative assembly method [29]. We use a com¬ 
mon glassy polymer as a model material: poly(methyl 
methacrylate) (PMMA, 120k g/mol) with fluorescent 
dye (Coumarin 153) for imaging. The ribbons are re¬ 
leased into a pool of water, at which point they sponta¬ 
neously form helices through a balance of surface tension 
and elasticity of the asymmetric cross-sectional geome¬ 
try. Details on helix formation and fabrication can be 
found in prior publications |28II30] . A micromanipulator 
with a carbon fiber attached at its end is subsequently 
used to fix one end of a helix and position it inside a 
microchannel (Fig. &) at the vertical center (ie. the 
center of the Poiseuille flow), which is 600 /im wide and 
100 /rm tall (Fig. Bi)- The carbon fiber is brought into 
contact with the helix and is fixed due to non-specific 
interactions. The flow velocity (v) is taken to be Vmax, 
where Vmax = Svavg/‘2 and Vavg = Q/A, where A is the 
channel cross-sectional area. 

The designed setup has the advantage of measuring 
both flow rate and helix geometry simultaneously in real¬ 


time with a flow sensor and a fluorescence microscope. In 
Figs. and [J, we present a typical flow cycle experi¬ 
ment to demonstrate helix shape recovery and flow con¬ 
trol. When the flow is turned on, the helix deforms along 
its helical axis in the direction of applied flow and in the 
absence of flow, the helix returns nearly to its original 
state. This particular helix is cycled three times from 
0 to 5 /rL/min (corresponding to u ~ 2 mm/s) and the 
flow rate history and the stretch ratio, A = H/Hq, are 
plotted along with corresponding micrographs. Here, Hq 
is the axial length of the helix in the absence of flow. At 
point b, A ~ 3.4 and recovers to point c where A ~ 1.1 
when the flow is turned off for 90 s. On the second cycle 
(point d), A « 3.4 and recovers to A « 1.15 and responds 
similarly in the third cycle, showing reversibility in our 
helices. The small, irreversible deformations observed are 
likely associated with creep deformations within the rib¬ 
bon material, but as shown below, these slight changes 
can be considered negligible for the focus of this work. 

Deformation of a helix in an external flow is caused 
by the hydrodynamic drag forces acting at each point 
along its contour length. Following resistive force the¬ 
ory, the drag force per unit length is given by: f = 
~C-L [v — (t • v) t] — C|| (t • v) t, where t is the local tan¬ 
gent of the ribbon backbone, v is the velocity of the 
fluid relative to the ribbon, and C_l and C|| are the fric¬ 
tional coefficients that define the resistance to motion of 
the surrounding fluid in the normal and tangential direc¬ 
tions from the ribbon, respectively [21E|. These frictional 
coefficients are proportional to the viscosity {rj) and a 
logarithmic correction dependent on the helical geome- 





































3 


try |55]. In general, for very elongated objects, the ratio 
C-l/C|| ~ 2. For the case of axial extension under flow, 
Kim and Powers [7] give an expression for the helix ex¬ 
tension in the limit that R/L is small and C-l/C|| ~ 2: 

( 1 ) 


where the velocity v is in the direction of the helical axis 
and the axial extension is defined as A7J = H — Hq. 

Guided by Eq. we measure AH as a function of v. 
From Fig. |^, the axial extension of the helices is non¬ 
linear with increasing flow velocity. We describe this 
non-linearity phenomenologically with non-linear finite 
extensibility [3T1132] which leads to: 



where the maximum extension is taken to be AH^ax = 
L — Hq. a typical experiment is shown in Fig. Helices 
with a range of sizes were created to examine the effects 
of helix and ribbon geometry; these range between R « 
3 — 15/im, p « 4 — 20/im, and L « 55 — 420/im. R is 
controlled by the bending stiffness B = El [55], where 
E = 2 GPa is a typical value for PMMA [33j. R and p 
are measured directly from the microscope images at zero 
flow rate and L is determined by the helical relationship 
L = -|-p2, where N is the number of turns. 

As expected, the assortment of helices display different 
flow-extension curves due to their varying shape and size, 
demonstrated in Fig. The dashed lines represent the 
fit of Eq. [^and our data are well fit to this relation. Using 
the corresponding R and L values at zero flow rate, Eq. 
leads to a best fit for B /C|| for all helices. Scaling with the 
measured helical geometries and determined B/C|| leads 
to a collapse of the data (Fig. [1)> validating the use 
of Eq. [ 2 I We determine the relationship between i?/C|| 
and i? ror flow experiments by fitting to the expression 
H/^ll = Cexp (aR), giving C = (1.7±0.7)xl0“^® m'^ s“^ 
and a = (3.5 ± 0.4)xl0® m“^ (Eig. [^). 

To determine frictional coefficients, we quantify B in¬ 
dependently with a recently developed micromechanical 
tool to measure the end-loaded force-extension relation¬ 
ship (Fig. [^) of our helices for different helix geometries 
(see refs. p^Si [34] for experimental details). Under end¬ 
loading conditions in the linear limit, the helical exten¬ 
sion is given by |7] : 


AH _ ER^ 


(3) 


At high extension, the force-extension relationship is non¬ 
linear and follows expressions developed previously by 
Pham et al. [31]|3S]. A plot of B vs. R provides the 
empirical relation B = Bq exp {f3R) with Bq = (2.6 ± 
l.l)xlO“^^ N m^, which is comparable to bacterial flag¬ 
ella, and /3 = (3.5±0.5)xl0® m“^ (Fig. [^). Importantly, 
we find a « /I, demonstrating that Cy is independent of R 




FIG. 2. (a) Fluorescent images of a helix with increasing 

flow velocities. At higher velocity, the helix begins to lose 
turns by rotating its free end. (b) Flow-extension curves for 
several helices, showing different extension due to the different 
helical dimensions {i.e. R and L). The second magenta data 
corresponds to (a). Dashed lines are a fit for a helix with 
non-linear, finite extensibility given in Eq. [^ (c) Data from 
(b) scaled by the helical geometries and determined B/C|| ■ (d) 
Semilog plot of B/C|| determined by the flow experiments as 
a function of the helix radius R. 
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FIG. 3. (a) Force-extension curve for an end-loaded helix, (b) 
The determined relationship between B and R in end-loaded 
experiments. 


within our experimental resolution and parameter range. 
Accordingly, a frictional drag coefficient can be quanti¬ 
tatively determined as C|| = -Bo/C = 1-5 ± 0.6 mPa.s. 

While different expressions of Cn have been proposed 
by different researchers mimiiiB], the general relevant 
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form for a circular cross-section is given by [26) : 


where q is usually taken as 0.09p and a as the radius of the 
cylindrical fiber itself. Since our helices’ cross-sections 
are not circular, but rather a shallow triangular ribbon 
(Fig. [^), we took the ribbon thickness to be the rele¬ 
vant length scale a. Although average ribbon dimensions 
are measured before transformation into helices, deter¬ 
mining the nanoscale cross-sectional thickness of specific 
ribbons in their helical form is not possible with the op¬ 
tical microscope used to record the helix deformations. 
Thus, the thickness t is determined for specific helices by 
relating R, measured optically, to established relations 
for B and measured aspect ratios of t/w by AFM and 
optical profilometry (Fig. Taking these values for t 
and the typical viscosity of water rj = 1 mPa.s, we calcu¬ 
late a theoretical C|| for each helix using Eq. providing 
C|| = 1.6 — 2.5 mPa.s. This range is in reasonable agree¬ 
ment with our experimental results: ('ll = 1.5±0.6 mPa.s. 
It must be noted that since the frictional coefficient (|| 
depends on the geometry of the helix as a logarithmic 
correction, the range of pitch and radius studied here are 
unfortunately not sufficient to resolve differences within 
our experimental resolution. 

Aside from their global extension, flexible helices dis¬ 
play non-uniform shape distributions when deformed in 
fluid flow, which has also been observed in helical flagella 
[Mils]. More specifically, it is observed that the turns 
are most stretched at the fixed end and continuously be¬ 
come less stretched along the helix approaching the free 
end. This is clearly visualized in an experiment of a long 
helix with several turns as shown in Fig. where L « 
320 fim and R « 4.5 ^m, and at an applied flow veloc¬ 
ity of u = 0.625 mm/s. Such non-uniform shapes are 
readily explained by distributed loading of helices: un¬ 
der flow, the force applied to a small element of a helix 
consists of local hydrodynamic drag on the element and 
the force accumulating along the helix from the free end. 
Mechanical equilibrium is then ensured by the equal and 
opposite force applied to the element by the rest of the 
helix that is further away from the free end. If we as¬ 
sume that the local hydrodynamic drag is independent 
of the position along the helix, as in the resistive-force 
theory discussed above, the total force applied to an el¬ 
ement of the helix from the free end grows linearly with 
the contour length s, measured from the free end. Con¬ 
sequently, the local pitch of the helix also grows linearly 
with s, as shown below. This situation is analogous to 
stretching of low-stiffness springs under gravity [SIEZI- 
Similarly, variations in the radius are observed when the 
local stretch is sufficiently high, as demonstrated near the 
fixed end of the helix in Fig. 

To quantitatively examine the shape distribution of 
our helices, we measure the local pitch p(s) by calculating 
the distance between the outmost points of neighboring 
turns along the helix. At zero flow rate, the pitch po is 




contour position, L-s (m) 


FIG. 4. (a) Helical ribbon with R ~ 4.5 /cm and L « 320 /rm 
in the absence of flow (top) and at n = 0.625 mm/s (bottom), 
(b) Pitch as a function of position on the contour length cor¬ 
responding to the helix in (a) for different flow velocity. The 
dotted lines are calculated from Eq. [^and the thin black rep¬ 
resents the calculated po- The error on pitch measurements 
is within the size of the points on the graph (< 1 pm). 


constant along the helix within small experimental varia¬ 
tions. Under flow, p{s) can be estimated as the difference 
of the axial displacement of the points s + l/2 and s — lj2 
obtained from Eq. 


p{s) = AH{s + 1/2) — AH{s — 1/2) + Po 
_ 2CiivR^ 

~ W 


-Is+Po, 


(5) 


where I is the contour length of one pitch, which we 
assume to be constant. This assumption is justified as 
long as no turns are lost during the experiment and irre¬ 
versible deformations are negligible. In Eig. |^, we plot 
Eq.j^with the corresponding values of (y, R and B deter¬ 
mined by our flow experiment for the different flow veloc¬ 
ities. In Eig. 1^, one can observe that for small velocities, 
the pitch vs. position dependence is well described within 
the experimental errors by Eq. Consequently, the ap¬ 
proximation of linear axial extension (Eq. holds at low 
velocities; however, deviations are observed at high flow 
velocities. Here the pitch distribution in the helix sec¬ 
tion experiencing the highest forces (closer to the fixed 
end) deviates from the linear approximation of Eq. In 
this part of the helix, measured pitches are smaller than 
predicted, corresponding to geometric strain stiffening at 
large extension also seen in Fig. [2>- Notably, our re¬ 
sults demonstrate a spatial manifestation of the crossover 
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FIG. 5. Qualitative example of a prominent helical transition 
from a global helical geometry to local uncoiled and coiled 
configurations at high velocity (of order ~ 10 mm/s). 


between linear and non-linear behavior of a helix under 
distributed loading. 

Finally, we interestingly observe large shape changes 
in helices at high velocities. Under these stronger forces, 
the helical shape uncoils to lose turns near the point of 
attachment (Fig. [2}t). More evident helical instabilities 
are observed through localized transitions of coiled to un¬ 
coiled helical geometry, as shown in Fig. These drastic 
deformations are reversible; when the flow is turned off, 
the helix relaxes to a shape that again is nearly identi¬ 
cal to the initial helix. Similar transitions of helical and 
straightened geometries have been observed in torque- 
free, end-loaded experiments on cholesterol helical rib¬ 
bons [38j as well as in rods of preferred curvature under 
gravity m- Moreover, similar conformational transi¬ 
tions at high flow rates are observed with flexible polymer 
chains [3S]. Thus, future studies will focus on these heli¬ 
cal transitions in viscous flow, which may provide impor¬ 
tant insight into unstable transitions that exist in helical 
systems found in nature. 


Overall, we have introduced a microscopic model sys¬ 
tem to measure the deformation, shape and frictional 
properties of flexible helices in low Reynolds number vis¬ 
cous flow and find that the global axial deformation is 
consistent with existing theory [7] . We demonstrate that 
with known ribbon properties and helical configurations, 
the shape distribution can be quantitatively predicted. 
Moreover, our experimental platform presents opportu¬ 
nities for theoretical advances on flexible helices in low 
Reynolds number flow; in particular, the effects of fluid 
viscosity or viscoelasticity, the friction and flow around 
deformable helices, the global-to-local helical shape tran¬ 
sitions, and potentially the effects of cross-sectional ge¬ 
ometry. Understanding these general helical behaviors 
both experimentally and theoretically will lead to funda¬ 
mental insights on natural helices, like flagella, as well as 
the development of synthetic helices, like swimming mi- 
crobots in fluid environments. Well-characterized helices 
can also be used to measure local forces in flows of sim¬ 
ple or complex fluids where local velocities can readily 
be measured by various techniques, like PIV or particle¬ 
tracking, while measuring local stresses presents a signif¬ 
icant technological challenge m- 
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